function H = Hessian(x)
    H=zeros(4,4);
    H(1,1) = 2 + 120 * (x(1) - x(4))^2 ;
    H(1,2) = 20 ;
    H(1,3) = 0 ;
    H(1,4) = -120 * (x(1) - x(4))^2 ;
    H(2,1) = 20;
    H(2,2) = 200 + 12 * (x(2) - 2 * x(3))^2 ;
    H(2,3) = -24 * (x(2) - 2 * x(3))^2 ;
    H(2,4) = 0 ;
    H(3,1) = 0 ;
    H(3,2) = -24 * (x(2) - 2 * x(3))^2 ;
    H(3,3) = 10 + 48 * (x(2) - 2 * x(3))^2 ;
    H(3,4) = -10 ;
    H(4,1) = -120 * (x(1) - x(4))^2 ;
    H(4,2) = 0 ;
    H(4,3) = -10 ;
    H(4,4) = 10 + 120 * (x(1) - x(4))^2 ;
end